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Abstract. We discuss various bifurcation problems in which two isolated 
periodic orbits exchange periodic "bridge" orbit (s) between two successive 
bifurcations. We propose normal forms which locally describe the corresponding 
fixed point scenarios on the Poincare surface of section. Uniform approximations 
for the density of states for an integrable Hamiltonian system with two degrees 
of freedom are derived and successfully reproduce the numerical quantum- 
mechanical results. 
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1. Introduction 

The periodic orbit theory [H [H [3] has made important contributions to the under- 
standing of quantum chaos [11 [5l [6] and to the semiclassical interpretation of quantum 
sheh effects in finite fermion systems [7|. Through semiclassical trace formulae, it 
relates the density of states of a quantum Hamiltonian to the sum over all periodic 
orbits of the corresponding classical Hamiltonian system. Gutzwiller's trace formula 
[1] assumes the periodic orbits to be isolated and therefore applies most directly to 
chaotic systems; it can, however, also be used for integrable systems with isolated 
orbits (see e.g., Oil]). 

In the derivation of the trace formula [1] , the stationary phase approximation is 
used for evaluating some of the trace integrals over the semiclassical amplitude, lead- 
ing to Gauss-Fresnel integrals. In systems with regular or mixed classical dynamics, 
periodic orbits can undergo bifurcations at critical values of a system parameter (e.g., 
energy or deformation). At such critical points, one or more of the Gauss-Fresnel 
integrals become singular and cause the divergence of the Gutzwiller amplitudes of 
the bifurcating orbits. This situation can be remedied [lOj by going to higher than 
second-order terms in the expansion of the action function, which appears in the phase 
of the trace integrand, around the critical point in phase space. The minimum number 
of terms required in this expansion are given by the so-called normal forms, which are 
characteristic for each type of bifurcation (see also jTl]) and lead to usually well- 
known catastrophe integrals. At the bifurcation points, one obtains in this way local 
uniform approximations for the semiclassical amplitudes that are finite and contain 
the contributions of all orbits involved in the bifurcation. 



Normal forms and uniform approximations for bridge orbit bifurcations 



2 



The local uniform approximations do, however, not reproduce the correct asymp- 
totic Gutzwiller amplitudes far away from the bifurcation point where all involved 
orbits are isolated. To achieve this goal, global uniform approximations must be devel- 
oped which interpolate smoothly from the local behavior at the bifurcation to the 
asymptotic regions of the isolated orbits. This has been done in Refs. [13l [Ml [15] for 
all generic bifurcations occurring in Hamiltonian systems with two degrees of freedom 
according to the classification of Meyer JJ^ (and listed also in [H]), in Refs. [T6l ITT] 
for codimension-two bifurcations, and in Ref. [18] for the transcritical bifurcation. (In 
passing, we mention that similar divergences of the Gutzwiller amplitudes occur when 
a symmetry is broken - or restored - under the variation of a system parameter. Local 
uniform approximations for symmetry breaking have been developed in Refs. [10[ I19j ; 
the prototype of a global uniform approximation for the breaking of U(l) symmetry 
was developed in Ref. [20], which inspired those mentioned above for bifurcations as 
well as global uniform approximations for the breaking of other symmetries [?!].) 

In this paper we investigate a type of bifurcation that has not yet been studied 
in this context and that we term bridge orbit bifurcation. Typically, this bifurcation 
consists of a pair of isolated orbits which are connected through a "bridge orbit" 
that only exists in a finite interval of the system parameter. Under its monotonous 
variation, the bridge orbit is born at a bifurcation of the first isolated orbit and then 
absorbed at a bifurcation of the second isolated orbit. This scenario has been found 
in both integrable and non-integrable Hamiltonian svstems[22l [23l |24] . It occurs, 
e.g., in the two-dimensional rationally deformed harmonic oscillator under a generic 
class of perturbations; other examples will be given in section [2l In the integrable 
case the bridge orbit forms a continuously degenerate family (i.e., a rational torus), 
while in non-integrable cases it is typically isolated. Since the two isolated orbits 
which exchange the bridge orbit typically are well separated in phase space, the 
bridge bifurcation is accompanied by global changes of the phase space structure which 
cannot be treated with the usual perturbative normal forms derived from the Birkoff- 
Gustavson expansion. This is different from the generic bifurcations of codimension 
one [TT] , and also from those of codimension two considered in [121 [H] j where only 
local changes in the phase space structure occur around a central periodic orbit that 
exists at all values of the relevant system parameter. 

In all the bifurcation types investigated so far in connection with uniform 
approximations, all orbits participating in the bifurcations become asymptotically 
isolated far enough from the bifurcation point (s). Even in the codimension-two 
bifurcations considered in [Tni [T7], where some isolated orbits only exist between 
two adjacent bifurcations, these have their well-defined partners which become 
asymptotically isolated at least on one side of one of the bifurcations. For the bridge 
orbits studied here, this is not the case. The bridge orbits are entirely intrinsic to 
the bifurcation scenario and cannot be linked to any 'external' orbit existing far away 
from the bifurcation. This causes a generic problem in the construction of the global 
uniform approximations which will be discussed in section [4] 

In section [2l we present some examples of bridge orbit bifurcations in integrable 
and non-integrable systems with two degrees of freedom. In section [3l we propose 
a new type of normal form (derived in detail in the appendix) which successfully 
describes the scenario of the bridge orbit bifurcation in integrable systems. In 
section [31 we derive local and global uniform approximations from this normal form, 
give analytic trace formulae for the semiclassical density of states and compare 
its numerically computed results with fully quantum-mechanical results for the 
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corresponding quantum Hamiltonians. Section [5] is devoted to a summary and 
concluding remarks. 

In the following sections, numerical integrations of the classical equations of 
motion are performed with Adams' Method (a kind of predictor-corrector method) to 
get solutions of high accuracy. Periodic orbits are obtained by searching fixed points 
(g,p) of the Poincare map in a suitable surface of section using a two-dimensional 
Newton-Raphson iteration method, whereby the stability matrix (linearized Poincare 
map) is obtained at the same time. The fixed points are smooth functions of the 
system parameter and can easily be followed through the bifurcation points under 
variation of the parameter, as well as new branches of fixed points emerging from the 
bifurcation points. 

2. Bridge orbit bifurcations in two dimensional potential models 

Let us start from a system with two degrees of freedom described by the Hamiltonian 

iJo(p,r) = ip^ + i|rr. (2.1) 

Since the potential is a homogeneous function of the coordinates r = {x,y), the 
Hamiltonian has the scaling property 

i7o(c5p,cir) = ci/o(p,r), OO, (2.2) 

and the equations of motion (EOM) are invariant under the following scaling 
transformation: 

p^c2p, r^c°r, t^c° 2i, (2.3) 

while the energy transforms as — > cE. Therefore, the phase-space profile is 
independent of energy and one has the same set of periodic orbits at all energies E. 
The action integral along a periodic orbit (po) has the following energy dependence 

SpoiE)^<f p-dr=f^y "<f p-dr = h£Tpo, (2.4) 

Jpo(E) \^0/ Jpo{Eo) 

where the dimensionless scaled energy £ and the dimensionless scaled period Tpo of 
the orbit are defined by 

whereby the reference energy Eq can be chosen arbitrarily. 

We now consider two kinds of perturbations of the system (|2.ip . The first one 
is introduced by a magnetic field perpendicular to the {x, y) plane. The motion of a 
charged particle in the plane is described by the Hamiltonian 

H, = ]^p' + ]^r'' ~ nr^~^{xpy-yp^), (2.6) 

where the radial dependence of the perturbation (with r = |r|) is determined such 
that the scaling invariance persists for any finite k. The strength of the magnetic field 
is proportional to Krt~^. Using polar coordinates r = (r, 0), p = (pr,Pe), (|2.6p reads 
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Figure 1. Scaled periods r of periodic orbits for the Hamiltonian 1 12.61 1 as 
functions of k. Isolated circular orbits and families of bridge orbits are indicated 
by dashed and solid curves, respectively. Circles represent bifurcation points. 
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Figure 2. Some short bridge orbits in Hamiltonian l|2.6|l for a = 2.4 and for 

different values of k. Upper row shows symmetric (1,1) bridge, middle and lower 
rows show asymmetric (2,1) and (3,2) bridges, respectively. 
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This system is integrable since pe = A is a constant of motion. For any non-zero k, 
there exist two isolated circular periodic orbits C± with different radii r± and angular 
momenta A±, which are found as solutions of the equations 

Kff(r,A)-ir" + |^-«;Art-i=i?, (2.8) 

^(r, A) = - r"-i- _ - .A (- - l) r^-^ = . (2.9) 

Evidently they depend on the parameter k. The radii r± correspond simply to the 
minima of the effective potential (|2.8p . Small oscillations around the circular orbits 
C± have frequencies ij}± for their angular and ri± for their radial motions which are 
given by 



^± = ^-^rl-\ n± = J^^{ri,Ai). (2.10) 

They become periodic for those values k± for which the two frequencies are 
commensurate, i.e., when 

^(.±) = ±^^^ii^, (2.11) 
UJ± n± 

with positive integers and rt_. Precisely at the values k — k±, the orbits C± 
must undergo bifurcations, because the trace of their stability matrix M equals two: 
Tr Mc_|_(k±) = +2. The interesting phenomenon now is that the two bifurcations are 
connected by a bridge orbit B that is created (or absorbed) at the bifurcations. It is 
actually a degenerate family of orbits with TrMB(K) = +2 for all values k_ < k < k^. 
More precisely, a bridge orbit B emerges from a bifurcation of the n_-th repetition 
of the orbit C_ at k = k_ and is absorbed at a bifurcation of the n_|_-th repetition 
of the orbit C+ at k = k+. It can therefore be labeled by the repetition numbers 
as B(n+,n_). Figure [1] shows the scaled periods of the shortest periodic orbits of 
the system (|2.6p as functions of the parameter k. At each crossing point of some 
repetitions of the two circular orbits, indicated by the pair of numbers (n-|_,7i_), one 
finds a bridge orbit family connecting them. The shapes of some of the bridge orbits 
are shown in figured) In figurel^l the scaled periods and the traces of stability matrices 
are plotted as functions of k for the circular orbits C± and the bridge orbits connecting 
them. Tt Mc_^{k) touches the horizontal line TrM = 2 at k = k±, and in between 
there exist the bridge orbit families having TrMe = 2. 

The second perturbation is introduced by an elliptic deformation. We modify the 
Hamiltonian (|2.ip as follows: 

Hp = lp' + lirMO)r, (2.12) 



fp{d) = Jijcos^e+-smH , 13 = ^i^^il . (2.13) 
V V ri+ 1 

The scaling rule (|2.2p persists for any /3, but the system is nonintegrable for (3^0. 
We take (3 as the deformation parameter, which is related to the axis ratio rj by the 
second equation in (|2.13p . For any /3 7^ 0, there are two isolated periodic orbits A^, 
and A^: straight-line librations along the x and y axis, respectively. Figure 0] shows 
the scaled periods of the shortest periodic orbits in the system (|2.12l) as functions of 
(3. Again one, finds bridge orbits B which bifurcate from the repetitions {n.j;,ny) of 
the orbits A^, and Ay near each crossing point. 
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Figure 3. Scaled period r (upper panels) and trace of stability matrix M (lower 
panels) of the symmetric (1,1) (left-hand side) and asymmetric (2,1) (right-hand 
side) bridge orbits for the Hamiltonian 1 12. 6| I with a = 2.2, shown by solid lines as 
functions of the parameter k. The dashed lines show the corresponding quantities 
for the circular orbits C±. 



In figure El we plot some of the shortest bridge orbits in the (a;,y) plane. The 
system (I2.12[) is non-integrable and the bridge orbits here are isolated. The symmetric 
bridges B(m, m) are bounded by two isochronous pitchfork bifurcations; the example 
of B(l,l) is illustrated in figure[6l In the left panel of figure[7l the graph of Tr M(/3) for 
the 2nd repetition of the A^; orbit touches the horizontal line Tr M = +2 at k = 0.628 
and two stable and unstable branches emerge from an island-chain bifurcation. The 
unstable and stable branches submerge sequentially into the Ay orbit at k = 0.688 and 
0.728, respectively, in two successive pitchfork bifurcations. For asymmetric crossings 
{ux, Uy) with Ux riy, the bridges consist of two isolated branches, i.e., one stable and 
one unstable orbit, which we label by Ti{nx, ny)s and B(nj;, ny)„, respectively. For the 
bridge orbits B(2,l), the common left end is a non-generic island-chain bifurcation of 
the 2nd repetition of the A^; orbit, while the right ends are two isochronous pitchfork 
bifurcations of the Aj, orbit occurring at two different deformations (see the left panel 
of figure [7]). The asymmetric B(n+,n_) bridges with n+,n^ > 2 also consist of 
two isolated branches with different actions and stabilities, but their end points are 
common at both ends of both branches, since they occur both at period- multiplying 
island-chain bifurcations (cf. the right panel in figure[7]for the B(3,2) bridges). 

In the limit a — s- 2, becomes the cranked and Hfj is the anisotropic harmonic 
oscillator. (Note that both systems have identical spectra after a suitable transfor- 
mation from a to k, cf. section 3.2.8 of [?])• In this limit, each bifurcation pair 
coalesces and the connecting bridge orbit shrinks to a point at which the two isolated 
orbits intersect. At the crossing points, one has locally periodic-orbit families of two- 
fold degeneracy due to the dynamical SU(2) symmetry of the rationally deformed 
harmonic oscillator in two dimensions. 

A scenario involving two bridge orbits is obtained if one breaks the U(l) sym- 
metry that the Hamiltonian (|2.13[) possesses at /? = 0. Let us e.g., modify the shape 
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Figure 4. Same as figure [T] but for the Hamiltonian 1 12.121 as functions of 
deformation parameter /3. Bridge orbits are isolated in contrast to those for (12.61 . 
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Figure 5. Some short bridge orbits in the Hamiltonian (|2.6|l for a = 2.4 and for 
several values of deformation parameters /3. The upper row shows the symmetric 
bridge B(l,l), the middle and lower rows display the stable and unstable branches 
of the asymmetric B(2,l) bridge. 
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Figure 6. Scaled period r (upper panel) and trace of the stability matrix M 
(lower panel) of the symmetric (1,1) bridge orbit for the Hamiltonian l|2.12| l with 
a = 2.2, shown by solid lines as functions of the deformation parameter /3. The 
dashed lines show the corresponding quantities for the Aj, and Ay orbits. 




Figure 7. Same as figure[6]but for asymmetric {nx,ny) bridge orbit bifurcations, 
with two pairs of values for the n^^-th repetition of the Ax orbit the and the riy-th 
repetition of the Ay orbit. The suffixes s and u represent the stable and unstable 
branches, respectively. 
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Figure 8. Same as figure |6] but for the shape function 112.141 1 with = ±0.002 
and a = 2.2. 



function fi3{9) in (|2.13|) by the following one: 

/;3,/34 {0) ^ cos2 6* + i sin^ 9- 13^ cos 46* , (2.14) 

with nonzero fi^. In this case, a second bridge orbit appears around f3 — for the 
symmetric {m,m) bifurcations. Figure [8] shows the properties of these bridge orbits. 
For /34 = it corresponds to figure [S] where the two isolated orbits and Aj^ intersect 
in a point with U(l) symmetry (/3 = 0) and their stability traces intersect at Tr M = 2 
(see the lower part of the figure). For f]^ ^ 0, this crossing point is split, so that the 
second bifurcation of each orbit occurs at different points on the (3 axis; the new pair 
of bifurcations is now connected by a second bridge orbit B2. For f]^ > 0, as shown on 
the left side of frgureO both bridges are stable, while for /34 < (see right side) one of 
them is stable and the other is unstable. All bifurcations here are of the non-generic 
pitchfork type. 

Finally, we mention a bifurcation scenario that has been discussed in |24j . Hereby 
a pair of isolated orbits exchange their stability via an isolated bridge orbit. Examples 
for this are found in the coupled quartic oscillator 

Ha = lp^ + \{x^ + y^) + axV. (2.15) 

Figure [9] shows a narrow region of the chaoticity parameter a. The shapes of two 
crossing isolated orbits F and P in the {x, y) plane are shown by inserts, as well as 
various shapes of the bridge orbit Q interpolating between those two shapes. Note the 
extremely small scale: the maximum value of | Tr M — 2| of the bridge orbit is smaller 
than 10~^. On a larger scale, the bridge orbit may not be observed numerically and 
the two isolated orbits F and P would appear to cross in a point. The bifurcation 
diagram then would look similar to that of a transcritical bifurcation [18] . The orbits 
F and P are created at a = 0.6315 in a period-tripling bifurcation from a straight-line 
libration along the y axis; at a = they become members of an integrable 3:2 torus 



Normal forms and uniform approximations for bridge orbit bifurcations 10 



2.e-07 



l.e-07 
-l.e-07 



-2.e-07 

0.005547872 0.00554805 0.005548228 

a 

Figure 9. Stability exchange of isolated orbits F and P in the coupled quartic 
oscillator l|2.15| l via two non-generic pitchfork bifurcations connected by an 
isolated bridge orbit Q. The inserts exhibit the shapes of the orbits in the (a;, y) 
plane. 

with U(l) symmetry. We refer to for details of this bifmcation scenario and to 
for details of the potential (|2.15p . A bridge bifurcation of the same type has 
been also found to occur in a two-dimensional spin-boson Hamiltonian |26| . 

3. Normal forms for some bridge-orbit bifurcation scenarios 

Normal forms are frequently used in singularity theory (see e.g., [27j ) and catastrophe 
theory (see e.g., [28]) to classify bifurcations. The natural variables of the normal 
forms for bifurcations in Hamiltonian systems with two degrees of freedom are the two 
canonical variables (g,p) spanning a projected Poincare surface of section transverse 
to the bifurcating parent orbit, and a bifurcation parameter e. These variables are 
usually chosen such that the bifurcation of the parent orbit occurs at (g,_p, e) = (0, 0, 0). 
The normal form function S{q,p,e) must fulfill the condition that its critical points 
correspond to the fixed points in (g,p, e) space in the neighborhood of the origin, and 
hence to the periodic orbits taking part in the bifurcation. This condition, however, 
is not sufficient to specify the normal form for a given bifurcation uniquely. Usually, 
S(q,p, e) is chosen to contain a minimum number of parameters and simple functions 
of the variables {q,p, e) - often just polynomial expressions in g, p and e - that yield 
the desired fixed-point scenario of a given bifurcation. For the generic bifurcations 
in two-dimensional symplectic maps according to the classification of Meyer [12] , the 
standard normal forms have been given in [10|, lllj . They can also be used for non- 
generic bifurcations of the same type (i.e., with the same fixed-point scenario in phase 
space) occurring in systems with discrete symmetries (see e.g., [29j). Normal forms for 
non-generic bifurcations in Hamiltonian systems with different fixed-point scenarios 
can be found in pE] (for the transcritical bifurcation) and in [16l [17] (for codimension- 
two bifurcations). 

For bifurcations involving more than two orbits, it is often useful to transform 




Normal forms and uniform approximations for bridge orbit bifurcations 11 

a=2A,K=0 







X 



Figure 10. Poincare surface of section for 112.61 1 with a = 2.4 and k = 0. 



the Poincare variables {q,p) to action-angle variables {ip,I) (cf. [TT| I15j): 



p = V2/cosv?, g = V2/sin(^, />0, ip e [0,2Tr) . (3.1) 

This becomes particularly useful when one considers integrable systems for which 
the normal form function does not depend on the angle ip, so that one only has 
to deal with a function S{I,e) depending on two variables. For instance, in [17], 
a simple bifurcation of a torus from an isolated orbit in the integrable Hamiltonian 
H = {p^ + Py)/2 + (x'^ + y^)/2 — could be described by the normal form 

r2 



Si{I,e) = Sa- el + ar 



(3.2) 



with suitably chosen constants 5*0 and a and bifurcation parameter e. The situation 
there corresponds to one half of that seen in figure [3] a torus B bifurcates from an 
isolated orbit C. The stationary condition for (|3.2[) is 

OS 'I , . 6 , , 

— (/B,.) = 0, ^ /B = -. (3.3) 

Since I must be positive definite, the torus B exists only for e/a > 0; its stability 
trace is always Tr Mb = +2 there. The isolated parent orbit arises from / = as a 
semiclassical end-point correction (see [T7] for details). Its stability trace is given by 
Tr Mc = 2 — e^, so that it is always stable except at the bifurcation point e = 0. 

In this section, we develop normal forms for some of the bridge-bifurcation 
scenarios described in the previous section. We start with the integrable model (|2.6p . 
Figure [To] shows a Poincare surface of section obtained for the parameters a = 2.4 and 
K = 0, fixing y = (and y > 0). Since the system is integrable, all fixed points lie on 
continuous curves which are intersections of rational tori with the {x,px) = iq,p) plane, 
except for the two isolated fixed points along the x axis which belong to the isolated 
rotating orbits C±. The corresponding bifurcation diagrams are seen in figure [1] where 
the rational tori correspond to the B(l,l) bridge orbit. The envelope of the curves in 
figure [TO] corresponds to the boundary of the classically accessible phase space given 
by Py = 0. 

The phase-space diagram seen in figure [TO] is completely analogous to that of a 
two-dimensional isotropic harmonic oscillator H = {p^. -f p^)/2 -I- ^^(a;^ -I- y^)/2. The 
tori correspond to the elliptic orbits, the isolated points to the circular orbits C± (with 
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both time orientations), and the hne a; = to the hbrating orbits. (The boundary 
contains the hbrating orbits along the y axis, i.e., with py — 0, which strictly cannot be 
seen in the Poincare plot.) It is useful to parameterize the orbits by their (conserved) 
angular momentum L = xpy — ypx- The circular orbits C+ and C_ have maximum and 
minimum value of L, respectively, and all intermediate nonzero L values correspond to 
elliptic orbits which are degenerate against rotations around the origin about an angle 
(?!> € [0, 7r). L = corresponds to the hbrating orbits which have the same degeneracy. 

This analogy suggests to use the angular momentum L and the rotation angle (j) 
of the harmonic-oscillator orbits to define action-angle variables ip = 2(f) and I — L/2 
and map these onto the Poincare variables of our present integrable system (|2.6p . 
This transformation is derived in [Appendix A[ The resulting mapping of {tp, I) to the 
Cartesian Poincare variables {q,p) is given by the relations (here for a; = 1) 

, — cos 9 

1 = V 2P y, ■ a = ' 
V 1 — sm cos if 



smOsmw „ „ 

2p-===^, p>0, 0<(p<27r, 
V 1 — sm u cos (f 

I = — q\J \p — — p^ — p cos , < 9 < tt , 



(3.4) 



which define an area-conserving canonical transformation {^p,I) ^ {(l,p)- Hereby 
is an angle parameterizing the angular momentum hy L — Lc cos 9, where Lc 
is the angular momentum of the C+ orbit in the harmonic oscillator. Note that 
the transformation defined by (|3.4p is more complicated than that given in (|3.ip : in 
particular, the action variable / here can have both signs. It is limited by the values 
±/9 yielding the fixed points {q,p) — (±-/2p, 0) which correspond to the two isolated 
orbits C± at their bifurcation points where the bridge orbit is created or absorbed. 

It turns out now that the bridge bifurcations in the integrable system (12.61) can be 
described by the same normal form function Si{I,e) as defined in p.2p above, except 
that here we have to use the definition of / in (|3.4p . As shown in [Appendix A[ p 
equals in the harmonic oscillator model. Here, p > is simply a parameter of 

the normal form that will be determined in section |31 

The stationary points of (/, e) are most conveniently found in terms of the 
variable 9. The stationary condition is 

— i =/9sin6l(e-2a/) = 0. (3.5) 
o9 

The stationary points satisfying sind — 0, i.e., = and tt, correspond to the two 
isolated circular orbits C± with / = ±p. In addition, there is another stationary point 
satisfying Ib — e/2a as in (j3.3p above, corresponding to the B torus. Since Ib must 
also fulfill the condition Ib = pcos9, it has real values only for 

- 2\a\p < e < 2\a\p , (3.6) 

so that the B torus only exists in the range of e values between the bifurcation points 
of the C± orbits. The actions of these periodic orbits are obtained by inserting their 
stationary values of / into the normal form (|3.2p : they become 

Sc± ^ Si{±p,e) ^ SoT<^P + ap'^, S'b = 5'i(/b, e) = - — . (3.7) 
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TrM = 



dpdq 



1 / d^S \ d^Sd'^S 
\ dpdq J dp^ dq^ 



S{q,P,e-) = Si{q,p,e) +qp, (3.8) 

whereby S{q^p,e) is the generating function of the Poincare map with initial 
momentum p and final coordinate q in the {q,p) plane at the value e of the bifurcation 
parameter. From this we obtain for our periodic orbits 

Tr Mc± = 2 - (e T 2ap)'^ , Tr Mb = 2 . (3.9) 

Hence we see that the normal form (|3.2I) with the mapping p.4p correctly describes 
the bifurcation scenario of the bridge orbits found in the model (|2.6p . as illustrated in 
figure [TJ 

In non-integrable systems, the normal form must depend also on the angle variable 
if. Let us consider the following normal form 

S2{I, if, e)^So-eI + al^ + b{p^ - I^) cos^ ip 

= S'o - epcos6' + p^{acos^ e + bsiT? Ocos^ (p) , (3.10) 

which respects the fact that the (^-dependent terms in any canonically invariant 
quantity (such as action, stability) should vanish for / = ±p. The stationary phase 
conditions become 
dS 

— — = - 2fep^ sin^ sin CO cos (i3 = , (3.11a) 

dtp 

dS 

^ psme(e-2pcose(a-bcos^ p)) =0. (3.11b) 

d9 

The above set of equations have the following solutions: 

sin6' = 0, 0±^O,TT ^ I±=±p, (3.12a) 

cos(p = 0, </5Bi = ^, ^ /bi = PC0s6Ibi = 7^ , (3.12b) 

2 2 2a 

sin<^ = 0, V3B2 = 0,7r Ib2 ^ pcos6>B2 ^ ^ ■ (3.12c) 

2{a — 0) 

The periodic orbits corresponding to the solutions (j3.12bp and (|3.12cp exist for 
— 2|a|p < e < 2\a\p and — 2|a — b\p < e < 2\a — b\p, respectively. The actions of 
these periodic orbits are 

S± = S'o T ep + o-P^, (3.13a) 

Sbi = S'o - — , (3.13b) 
4a 

^B2 = So- ^, ,, +bp^, (3.13c) 
4(a — o) 

and the traces of their stability matrices are 

TrM± =2-(e=F2ap)[e=F2(a-fe)/9], (3.14a) 

TrMBi =2 + -(e-2ap)(e + 2ap), (3.14b) 
a 

TrMB2 = 2 ^—[e-2{a-b)p][e + 2{a-b)p]. (3.14c) 
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Thus, the periodic orbits l|3.12b[) and (|3.12cp are the two bridge orbits which connect 
the two periodic orbits corresponding to (|3.12ap at the bifurcation points Cbif — ^^ap 
and ±2(a — b)p. The bridge orbits found in figure [5] can be regarded as of this type. 
For b — a (or a = 0) the bridge orbit B2 (or Bl) shrinks to a single point e = 0. 
This corresponds to the situation found for the (1,1) orbit in figure [6l The orbit Q in 
figure [9] can be regarded as the bridge orbit B2 in the hmit of extremely small |a — 

Figure [TT] shows the Poincare surface of section for the Hamiltonian (|2.12p , with 
deformations given by ()2.14p . for a = 2.2, P4 = 0.002 and several values of /3 along 
the right half of the (1,1) bridge bifurcation region (cf. the left panel of figure [5]). The 
Poincare variables {q, p) are here defined by g = (x + y)/ V2, p = (px+Py)/ at the 
surface of section x — y = 0, i — y > 0. At /? = 0, the existing periodic orbits are the 
bridge Bl at {q,p) ~ (0, ±1), the bridge B2 at {q,p) = (0,0), and the diameter orbits 
Aa;, Ay at {q,p) « (=f1,0). Both bridges are stable and the diameters are unstable. 
At the bifurcation point /3 = 0.016, the bridge B2 merges into the A^; orbit and A^; 
becomes stable. At the bifurcation point (3 = 0.077, the bridge Bl merges into the 
Ay orbit and Ay becomes stable. By fitting (I3.13P and (|3.14p to the periodic orbit 
quantities for a — 2.2 and /34 = 0.002, we obtain ap = —0.12 and bp = —0.145. 
The bifurcation points are thus e = ±0.24 and ±0.05 for the Bl and B2 bridges, 
respectively. Figure [12] shows the contour map of the normal form (|3.10p with the 
above values of a and b, and with e corresponding to each value of /? in figure [TT] One 
will see that the phase space profiles of figure [TT] are nicely reproduced by this normal 
form. 

For the asymmetric bridges (71+ , ) with n+ ^ n_ we have so far not found any 
suitable normal form. One may need a mapping based on the anisotropic (rational) 
harmonic oscillator instead of that given in (|3.4p which is based on the isotropic 
oscillator. This will be investigated in further research. 

4. Uniform approximations for the density of states 

4-.1. Local uniform approximation 

In this section we evaluate the semiclassical level density around the bifurcation points 
using the normal form obtained in the above section. In the following wc limit ourselves 
to the integrable model described by the Hamiltonian (|2.6p . The calculation of the 
parameters in a normal form from a given Hamiltonian is in general a very difficult 
problem. The idea of the uniform approximations is to avoid their direct calculation 
by relating them to the local invariant properties of the participating periodic orbits 
in the vicinity of a bifurcation, which can be obtained numerically from a periodic- 
orbit search. The bridge bifurcations occurring in the Hamiltonian ()2.6p involve three 
orbits: the two circular orbits C± (with repetition numbers n+ and n_) and the bridge 
orbits B(n+,n_) (cf. figures [T] [3]). From these, we can determine five independent 
quantities: the three actions S± and 5b, and the two stability traces Tr M± (recall that 
TrMe = 2 is constant). The normal form given by p.2p and (|3.4p contains the four 
parameters ^o, e, a and p, which we can determine using four of the above five orbit 
properties. In problems with only one bifurcation point, such a way of determining 
the normal form parameters does not contradict with remaining unused quantities, 
but in the present situation of the bridge bifurcations, the values of the parameters 
do depend on which quantities are used. This problem is related with the global 
nature of the bridge bifurcations, where the parameters undergo significant changes 
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Figure 11. Poincare surface of section for the Hamiltonian 1 12. 121 1 with 
deformation l|2.14ll . plotted for a = 2.2, /34 = 0.002 and several values of /3 in 
the (1,1) bifurcation region. (See text for definitions of the Poincare variables 




Figure 12. Contour map of the normal form Il3.10| l with ap = —0.12, bp = —0.145 
and several values of e in bifurcation region corresponding to each panels of 
figure [TT] (q,p) are plotted in units of y/2p. 
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between the two bifurcation points. Different from all bifurcations treated so far in 
the literature in uniform approximations, the bridge orbit here has no 'external link', 
i.e., there is no external orbit outside the bifurcation interval to which its properties 
can be asymptotically linked. This leads to a slight ambiguity in determining the 
parameter o, as we shall see below. 

In order to determine the normal form parameters uniquely, at least locally for 
each given e, we add one more term to Si{I,e) in (|3.2|) and use the normal form 

S•3(/,e) = 5o + e/ + a/2 + 6/^ I ^ psm9 . (4.1) 

The stationary-phase analysis then predicts the properties of the periodic orbits to be 



S± = So + ap^ T (ep - bp^) 



Sn — So — 



'■{1 + 2be/a^ + y/l + 3be/a^) 



(4.2a) 
(4.2b) 



a(l + v/l + 36e/a2)3 

TrM± = 2- [e- (±2ap + 36p^)]^ (4.2c) 

The five normal form parameters (5*0, e, a, b, and p) are now uniquely determined by 
the five equations in (|4.2|) . although these cannot be solved analytically. Due to the 
scaling rules, the parameters have the following energy dependences: 

b 



So = TlSto , p = hSp , 



(4.3) 



where tq, p, a, and b are dimensionless constants. Ideally, these four parameters should 
not depend on the bifurcation parameter e throughout the bifurcation region. 

Figure fT3l shows their results which we have determined numerically for the (1,1) 
and (2,1) bridge-orbit bifurcations for a — 2.02. In the center panels, we show besides 
To also the scaled periods r of all three periodic orbits. The parameters a, b, Sq 
and p turn out to be approximately constant throughout the bifurcation region, as 
hoped, in particular for the symmetric (1,1) bridge. Note that at e = 0, the value 
of b is exactly zero for the symmetric bridge bifurcation, so that the cubic term in 
(|4.1[) does not contribute there. Furthermore, the contribution of the cubic term bl^ 
to the actions of the periodic orbits remains much smaller than that of the quadratic 
term al^ throughout the whole bifurcation region in both cases (1,1) and (2,1). It will 
therefore be omitted again in the following. 

The combined contribution of all orbits involved in the bifurcation to the semi- 
classical level density is given by [9l [TO] 

^^Re Jdqjdp'i{q,p,e)exp ^ {5'(g,p, e) - gp} - yi/ , (4.4) 



dg{E) 



where S{q^p,e) is the generating function defined in (j3.8p which contains the appro- 
priate normal form, and the amplitude function ^!{q,p, e) is given by 

1/2 

(4.5) 



^ 9S 



d^S 



dqdp 



Here v is the Maslov index of the bridge orbit labeled by [n^ , ^^- ) • It is given by 

V — 2nr , rir — rt+ + 7i_ , (4-6) 

where Ur is equal to the number of librations in the radial direction. The amplitude 
function "^{q,p, e) usually has a moderate dependence on the variables q,p, e and can 
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Figure 13. Normal form parameters determined from periodic orbit quantities. 
Left and right panels are results for the (1,1) and (2,1) bridge bifurcations, 
respectively, for a = 2.02. In the center panels, the scaled periods r of all three 
orbits are shown besides tq = So/hS. 



be replaced by its value at the origin, ^'(0, 0, 0) = dSo/dE = To, which is the average 
period of the orbit cluster. By transforming the variables from {q,p) to (</?,/) and 
integrating over ip, one obtains 



Sg{E) ~^Re e'iSo/n-../2) f" ^^^n-ei+ar 



)/h 



(4.7) 



Here, we omitted the cubic term bP in the normal form (|4.ip as stated above. The 
integral on the right-hand side can be analytically expressed as 



J' = e-'/^'-'' {a+{c+ + is+) + a_(c_ + i.s_)} , (4.8) 



c± = C(|a;±|), s± = S(|a;±|), ^±^Y^(^T^)' <^±=sgnx±, (4.9) 

where C(a;) and S(x) are the Fresnel functions defined by 

S{x)^J sm[^t^jdt, C{x) = J cos[^t^jdt. (4.10) 

Inserting (|4.8p into ()4.7p . one obtains the local uniform approximation to the level 
density: 

Sg{E) = /° Re e'^^-/'^-''/^) {a+{c+ + is+) + a^{c^ + is^)} . (4.11) 
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In terms of the scaled dimensionless normal form parameters given in (j4.3p . the scaled 
level density becomes 

= \f^ToRee'^^-'-^''^'Ha+ic++is+)+a^{c^+is^)}, (4.12) 
V zTra 

where tb = Sb/Tt-^ — + /Aa is the scaled period of the bridge orbit, c± and s± 
are given by (|4.9p with the arguments 




(4.13) 



4-2. Global uniform approximation 

The local uniform approximation (|4.1ip is valid only in the vicinities of the bifurcation 
points, where the properties of periodic orbits are nicely described by the normal form 
(j3.2p . Far from the bifurcation points in the e variable, the standard asymptotic trace 
formulae should work, which are given in the form of the Gutzwiller formula [I] for the 
isolated orbits and the Berry- Tabor formula [30] for the family of bridge orbits. The 
purpose of the so-called global uniform approximations is to interpolate between the 
local result (14. lip and the standard trace formulae. For the generic bifurcations, such 
global uniform approximations were developed by Sieber and Schomerus [T31 HH US] . 
We shall presently follow their procedure to derive a global uniform approximation for 
the integrable bridge bifurcations. 

In order to go beyond the local uniform approximation, one can start from 
equation (|4.4p . but one has to expand the amplitude function '^{q,p,e) around the 
bifurcation point in a way similar to the normal form function S{q,p, e) for the action 
integral. We find that it is sufficient here to take into account only a linear term in /. 
We thus write 5g{E) as 



ITT 
— V 

2 



(4.14) 



^^^^f d/(« + PI) exp ^ [S{q,p) - qp] 

In principle, the parameters a and (3 are determined from higher-order expansion 
coefficients in the normal form (cf. |14p. but in practice one may determine them 
from the condition that the integral (|4.14p reproduces the asymptotic contributions of 
the bifurcating orbits to the standard trace formulae far from the bifurcation points. 
There, the action differences between different periodic orbits corresponding to the 
stationary points of the normal form p.2p are much larger than h. This is equivalent 
to taking the asymptotic expansions of the Fresnel integrals when their arguments are 
much larger than unity. 



, 2(ap2^£p + eV4a) 2 S±-Sb ^ ^ 

"^ = V -h = V^^^»'- 

Their asymptotic forms are (see e.g., |31j ) 

L(x) ~ — I sm —x , Six) ~ cos —x , 

^ ^ 2 nx 2 ^ ' 2 nx 2 

C(a;)-f iS(a;) ~ ^ e'"" x>l. (4.15) 

V2 TTX 
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1 



Re 



+ ^ /'gi(S+A-7r(i.+ l)/2) _ gi(S-A~7r(,.+ l)/2) 

2a \ 



-)} 



(4.16) 



with 5b and S± given by p.7p . Inserting (I4.15p . one has 



£/3 



7r?i yjixajh, 

7r?i|2ap — e| 
a — I3p 



cosl---.- 



cos I 77 + O--) 



(4.17) 



The first term on the right-hand side can be identified as the contribution of the bridge 
orbit with its Berry- Tabor amplitude 



A 



BT = — 2 
TTft 



JC 



dP 



and we have 



(4.18) 



(4.19) 



Tb being the period of the primitive bridge orbit. The second and third terms in the 
right-hand side of (14.17^ are identified as the contributions of two isolated orbits with 
their Gutzwiller amplitudes 

A± = (4.20) 

7rnn±i/|2-TrM±| 

where T± are the full periods of these orbits and n± are their repetition numbers. 
Furthermore one finds 

a±(3p^T±, 



2ap±e = a±n±y/\2 - TrM±| . 
From (|4.20p and (I4.2ip . we can determine the parameters a and /? by 



a-i(r++T_) 



2wha £±I± + 

A+ ' A- 



(4.21) 



(4.22) 



Inserting them into (|4.16p . we obtain the global uniform approximation for the level 
density 



{<7+{c+ + is+) + (7_(c_ + is-)} 



a+T+ a^Tl 



54 



COS -f - T7(^ + 1) - COS —--(;/ + 1) 



(4.23) 
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with the arguments of the Fresnel integrals given by 




(4.24) 



Note that all quantities entering this formula are determined from the invariant 
properties of the periodic orbits. 

The global uniform approximation (|4.23p is the main result of this section. We 
shall refer to it in the following as UA2. It becomes important in particular when 
the two bifurcations, at which the bridge orbits are created and annihilated, are close 
to each other so that neither of the isolated circular orbits are in the asymptotic 
region; the two bifurcations then cannot be treated separately. Sufficiently far from 
the bifurcation points, where all orbits reach their asymptotic domains, one can use 
the asymptotic trace formulae, referred to as ASY in the following, i.e., the Gutzwiller 
formula [T] for the isolated circular orbits and the Berry- Tabor formula [W for the 
bridge orbit families. If this is also the case for a central region between the two 
bifurcations, where the bridges are present, their bifurcations from/into the isolated 
orbits can be treated separately in the uniform approximation given in [17 , which we 
shall refer to as UAl. 

4-3. Numerical results 

In the following, we test numerically the various semiclassical approximations to the 
level density. We restrict ourselves to the coarse-grained level density in the scaled 
energy variable £, which quantum- mechanically is defined in terms of the scaled energy 
spectrum {£i\ by a sum over normalized Gaussians with width 7: 



In the semiclassical trace formulae, each contribution of a periodic orbit (po) then has 
to be multiphed by an exponential damping factor exp{ — (7rpo/2)^} (see e.g., [7]), so 
that the contributions of longer orbits are suppressed and the sum over the periodic 
orbits converges. In our numerical results given below, we have used 7 — 0.3 in all 
cases. With this value of 7, the condition for the above damping factor to be smaller 
than 10^^ corresponds to Tpo > 14, whose contributions can be safely neglected. 

In figure [141 the oscillating part of the scaled level densities dg{£), labeled by ASY 
and UAl, are compared with the quantum-mechanical result (QM) for the system (|2.6p 
with a = 2.4 and n = 0.31. These values of the parameters correspond to a point 
close to the left bifurcation of the bridge orbit B(2,l). In calculating the semiclassical 
level densities, we take into account the circular periodic orbits C"^ with repetition 
numbers < 3, the orbit C_ and the bridge family B(2,l) (see figure [1] for the 
values of their scaled periods). We see that the uniform approximation UAl for 
the bifurcating circular orbits with the bridge orbit improves the semiclassical level 
density over the asymptotic one and nicely reproduces the quantum results near the 
bifurcation points. 

Figure [TSl shows the result of the global uniform approximation UA2 for a = 2.02. 
K = 0.0025 and k = 0.322 correspond to the bifurcation points of the bridge orbits 
B(l,l) and B(2,l), respectively. In the periodic orbit sum, all periodic orbits with 
T < 15 are included; namely, C!^ and bridge B(n, n) with n < 3 for k = 0.0025, 
and C"^ with n+ < 3, C_ and B(2,l) for k = 0.322. In comparison to the separate 
uniform treatment of the two bifurcations (approximation UAl), the UA2 formula 




(4.25) 
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Figure 14. Oscillating part of the scaled level density Sg{£) for the Hamiltonian 
112. 6| I with a = 2.4 and ft = 0.31. Solid, dashed and dotted curves represent UAl, 
ASY (divided by 10) and QM, respectively. In the semiclassical level densities, 
the periodic orbit sum includes the orbits with n+ < 3, the orbit C_, and 

the bridge B(2,l). 




Figure 15. Comparison of uniform approximations for a = 2.02. Solid, dashed 
and dotted curves represent UA2, UAl and QM, respectively. Upper and lower 
panels are calculated at the bifurcation points of symmetric bridge (1,1) and 
asymmetric bridge (2,1), respectively. In the upper panel, the periodic orbit sum 
in the semiclassical level density is taken over C"* and the bridge B(n^_, n_) with 

n± < 3, while in the lower panel, C^"*" with n+ < 3, C— and the bridge B(2,l) 
are taken into account. 



(j4.17p reasonably improves the level density, but we obtain a slight overestimation. Let 
us consider the origin of this deviation. We note that for a ~ 2, there is no asymptotic 
region between the two bifurcation points. Therefore, the procedure to determine the 
parameters to reproduce the Berry- Tabor asymptotic form is not justified. In the 
UA2, the normal form parameter a is determined such that the asymptotic form 
reproduces the Berry- Tabor trace formula. This corresponds to using the curvature 
K, of the torus for the normal form parameter a as in (|4.19|) . Figure [16] compares 
the values of a determined by the two methods, i.e., by (|4.19p and (|4.2p . They are 



Normal forms and uniform approximations for bridge orbit bifurcations 



22 



uniform 

(X= I.KjZ. c II rvature 



(1,1) 



a=2.02 
(2,1) 



uniform 
curvature 



0.04 0.06 



Figure 16. Normal form parameter a = h£a (see 14.311 ) determined from 114.21 1 
(solid curve) and that from ||4.19| | (dashed curve). Left and right panels show the 
results for (1,1) and (2,1) bridge orbits, respectively. 



significantly different from each other, and the global uniform approximation (UA2) 
does not coincide with the local uniform approximation (|4.1ip in the bifurcation region. 

Using the local uniform approximation (|4.1ip with the normal form parameters 
given by (|4.2p (referred to as UA2L), the results are much improved near the 
bifurcation points as shown in figure fTTi but at the middle of the bridge, the quantum 
results lie between UA2 and UA2L. 

Our numerical results can be interpreted as follows. The normal form (j3.2p or 
(|4.ip is constructed to reproduce the bifurcation properties of the three participating 
periodic orbits. Therefore, it is most reliable around the bifurcation points e ~ Cbif. 
However, the bridge orbit undergoes large changes between the two bifurcation points 
which are globally separated in the phase space, and therefore higher-order terms in 
I in the normal form could contribute significantly in the middle region of the bridge. 
In fact, the curvature /C of the torus, which is an invariant property of the periodic 
orbit family and important for determining its contribution to the level density, is not 
correctly described by the form (|4.ip . as shown in figure [HI When the parameter a is 
shifted toward the value determined by the curvature /C, the agreement between the 
local uniform approximation and the quantum results becomes better. This could be 
achieved by normal forms with higher-order terms in / which, however, would render 
the global uniform approximation more complicated and less analytic. 

5. Summary 

We have investigated the appearance of bridge orbits, which connect two isolated 
orbits via two successive bifurcations near the points where their periods and stabilities 
coincide, in various Hamiltonian systems with two degrees of freedom. In a class of 
integrable systems, the bridge orbits form degenerate families. For these we used 
a mapping derived from the Poincare variables of the isotropic harmonic oscillator 
to derive a very simple normal form from which all the invariant properties of the 
participating periodic orbits can be derived analytically. Using this normal form, we 
have derived analytical uniform approximations for the semiclassical level density of 
the corresponding quantum systems. Although the normal form parameters could not 
be determined uniquely and their values undergo slight variations between the two 
bifurcations, the numerical agreement between the semiclassical and the quantum- 





Figure 18. Same as figure (TT] but calculated at the middle of two bifurcation 
points of (1,1) and (2,1) bridge orbits. 
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mechanical coarse-grained level densities is very satisfactory. 

We expect that the remaining differences and the slight variations of the normal 
form parameters can be reduced by including more terms in the normal form. 
This would be at the cost of losing the simple analytical forms of the uniform 
approximations and of having to determine more parameters numerically. 

The exploration of suitable normal forms for bridge-orbit bifurcations in non- 
integrable systems, such as those shown in figures [7] and [H is the subject of further 
studies. The uniform approximations for non-integrable system form also an important 
subject for understanding the deformed shell structure in realistic nuclear mean-field 
models [23 [23]. 
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Appendix A. Derivation of the mapping for bridge orbit bifurcations 

In this appendix, we derive a mapping from the Poincare variables (g, p) to action-angle 
variables (iy9, /) that are suitable for the normal forms of bridge orbit bifurcations. Let 
us consider the isotropic two-dimensional harmonic oscillator 

H=\{pI+pI) + \l.\x' + V^). (A.l) 

All its trajectories are periodic with period T = 2tt/lo] they are ellipses which may 
degenerate to a circle or to linear librations. We parameterize these periodic orbits 
using the two constants of motion energy E and angular momentum L = xpy — ypx ■ 



The orbit whose longer semiaxis lies on the x axis is written as 

xo{t) = qc [cos ^ + sin cos{ujt) , (A. 2) 

ya{t) = 5c f cos - - sin - j sin(cji) , (A.3) 

Qc V-^c/w, Lc = E/lj, L^Lc cos 9. (A. 4) 



Here Lc is the maximum angular momentum at fixed energy, which is that of the 
circular orbit running anti-clockwise in the {x,y) plane. 

Rotating this orbit by an angle ip/2 about the origin, one obtains the most general 
orbit 

x{t) = xo{t) COS ^ - yo(t) sin ^ , 

y{t) = xo{t)sm^ + yo{t)cos^, 0<(^<2^. (A.5) 

The projected Poincare surface of section E := {{x{ti),px{ti))\ y{ti) =0, y(ti)>0} at 
successive times [i — 1,2,...), with ti+i = ti+T, defines the Poincare variables 
(QtP) = {x{ii)-,Px{ti)). Each point corresponds to an orbit with period T, 

and therfore the Poincare variables {q,p) have a one-to-one correspondence with the 
variables {9,ip), given by 
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The Jacobian of the transformation {if, 9) — > {q^p) is 

We now define the action variable / by 

I = pcose, p^^pcqc = ^Lc, -p<I<p. (A. 8) 

Comparing to (|A.4[) . one sees that / = L/2. One also sees easily that the trans- 
formation {ip,I) — > {q,p) becomes canonical since 

d{q,p) 



d{ifj) 

For the inverse transformation, one finds from (IA.6P and (lA.SP the relation 



(A.9) 



^= ^9v/4pcc;-p2_^2q2^ ]^[pi+uj''q^)<2pu;{^E). (A.IO) 
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